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ABSTRACT 

We consider the hypothesis that galactic magnetic fields are primordial. We also 
discuss the various objections to this hypothesis. We assume that there was a magnetic 
field present in the galactic plasma before the galaxy formed. After the galactic disk 
formed, the lines of force thread through it and remain connected to the external 
cosmic medium. They enter through one side of the disk, proceed horizontally a 
distance I in the disk, and then leave through the other side. We find that the lines of 
force are stretched rotation of the galactic disk, which amplifies the toroidal component 
of the field and increases I. When the magnetic field is strong enough, it produces 
ambipolar velocities that try to lift the line out of the galactic disk but in opposite 
directions on different parts of the line. The result is, instead of the line being expelled 
from the disk, its horizontal length I is shortened, both in the radial, and in the 
toroidal direction. This leads to a reduction of the rate of horizontal stretching, and, 
finally, a reduction in the magnetic field strength. After a sufficient time, the magnetic 
field at all points goes through this stretching and reduction, and the field strength 
approaches a universal function of time. This function is slowly decreasing, and only 
depends on the ambipolar properties of the interstellar medium. At any given time 
the magnetic field is toroidal, and has the same strength almost everywhere. On the 
other hand, it turns out that its direction varies rapidly with radius, changing sign 
every one hundred parsecs field dominate over that of the other. The resulting field 
has a net Faraday rotation. If such a field were observed with low resolution in an 
external gal axy then the field would appear toroidal in between the spiral arms. The 
spiral density wave would turn it so that the lines appear to trace out the spiral arm, 
although the apparent lines really are the sum of pieces of magnetic lines as they cross 
the disk. They do not necessarily extend very far along the arms. We contend that 
this model of the magnetic field, which arises naturally from a primordial origin, can 
fit the observations as well as other models for the magnetic field, such as those arising 
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from the mean field dynamo theory. Finally, because the field lines are topologically 
threaded through the disk they cannot be expelled from the disk. This counters the 
objection against the primordial origin, that such a field could not survive very long in 
the galaxy. 

Subject headings: galaxies:magnetic fields-ISM:magnetic field— 
magnetohydrodynamic:MH-plasmas 
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1. Introduction 

There are two schools of thought as to the origin of the galactic magnetic field. The first 
school holds that, after formation of the galactic disk, there was a very weak seed field of order 
lCT 17 gauss or so, and that this field was amplified to its present strength by dynamo action 
driven by interstellar turbulence. The description of how this happens is well known (Ruzmaikin, 
Shukorov and Sokoloff 1988), and is summarized by the mean field dynamo theory. 

The second school of thought concerning the origin of the galactic field is, that it is primordial. 
That is, it is assumed that there was a reasonably strong magnetic field present in the plasma 
before it collapsed to form the galaxy. It is assumed that the field was coherent in direction and 
magnitude prior to galactic formation. The actual origin of the magnetic field is not specified. 
However, it could be formed by dynamo action during the actual formation of the galaxy, (Kulsrud 
et al 1996). 

This paper discusses the evolution of the structure and strength of such a primordial magnetic 
field during the life of the galaxy. The origin of the field prior to formation of the galaxy is not 
considered. 

It was Fermi(1949) who first proposed that the galactic magnetic field was of primordial 
origin. Piddington(1957, 1964, 1967, 1981) in a number of papers suggested how this might 
actually have happened. However, he supposed the magnetic field strong enough to influence the 
collapse. One of the first problems with the primordial origin is the wrap up problem. It was 
pointed out by Hoyle and Ireland(1960, 1961) and by Oke et. al. (1964), that the field would wrap 
up into a spiral similar to the spiral arms in only two or three galactic rotations. It was supposed 
that this is the natural shape of the magnetic field lines. If the winding up continued for the fifty 
or so rotations of the galaxy, the field lines would reverse in direction every one hundred parsecs. 
This seemed absurd. The various attempts to get around this problem involved magnetic fields 
which were strong enough to control the flow, and also strong radial outflows. Several forceful 
arguments against the primordial origin were advanced by Parker. The arguments were that the 
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field would be expelled from the galaxy, either by ambipolar diffusion (Parker 1968), or by rapid 
turbulent diffusion (Parker 1973a, 1973b) 

In this paper we reexamine the problem of a a primordial field. We proceed along lines first 
initiated by Piddington. Just as in his model for the primordial field, we start with a cosmic field 
with lines of force threading the galaxy. We further assume that, after the collapse to the disk, the 
lines remain threaded through the galactic disk. The lines enter the lower boundary of the disk 
in a vertical direction, extend a short horizontal distance in the disk, and then leave through the 
upper boundary. Thus, each line initially extends a finite horizontal distance in the disk. 

However, in contrast to Piddington's model, we assume the field is too weak to affect the 
plasma motions, especially the rotational motion about the galactic center. Consequently, it tends 
to be wrapped up by the differential rotation of the galaxy. In addition, we include additional 
physics of the interstellar medium in which the magnetic field evolves. After toroidal stretching 
strengthens the toroidal field sufficiently, a strong vertical force is exerted on the ionized part of 
the interstellar medium forcing it through the neutral part. We contend that this force will not 
expel the entire lines of force from the galactic disk, but the resulting ambipolar velocity will 
only tend to shorten the length which each line spends in the disk as it threads through it. In 
particular, ambipolar diffusion will decrease the radial component of the horizontal field and 
shorten its radial extent in the disk. This will decrease the toroidal stretching of the field line, 
and as a result the magnetic field strength will approach a slowly decreasing saturated condition. 
Thus, after several gigayears, the lines will end up extending a longer distance in the azimuthal 
direction than they did initially, but a much shorter distance in the radial direction. At any given 
time after the field saturates, the final strength will be independent of its initial value. In addition, 
the field strength will depend only on the ambipolar properties of the interstellar medium. [For 
earlier work see Howard (1995) and Kulsrud(1986, 1989, 1990).] 

As a consequence of ambipolar diffusion plus stretching, the magnetic field will be almost 
entirely toroidal, and it will have the same strength everywhere at a fixed time. However, its sign 
will depend on the sign of the initial radial component of the field at its initial position. Because 
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of differential rotation the toroidal field will still vary rapidly in sign over a radial distance of 
about one hundred parsecs. 

This model for the magnetic field evolution would seem at first to leave us with a wrap up 
problem and produce a field at variance with the observed field. However, the initial regions 
where the radial field is of one sign are expected to be of different area than those of the other 
sign provided that the initial field is not exactly uniform. As a consequence, at the present time 
in any given region of the galaxy, the toroidal field should have a larger extent of one sign than 
of the other, even, though the sign varies rapidly. If one now averages over larger regions than 
the size of variation, one would see a mean magnetic field of one sign. This averaging is actually 
performed by the finite resolution of the observations of the magnetic field in our galaxy or in 
external galaxies. From the observations one would not be aware of this rapid variation. 

If one ignored the spiral arms, the field would be almost completely azimuthal. However, it 
is known that the density compression in the spiral density wave, twists the magnetic field to be 
parallel to the spiral arm and increases its field strength (Manchester 1974, Roberts and Yuan 
1970). In observing radio emission from external galaxies, one tends to see the radiation mostly 
from the spiral arm, where cosmic rays are intense and the magnetic field is stronger. Because 
in these regions the field is aligned along the arms, one would naturally get the impression that 
the field extends along the entire arm. In fact, on the basis of our model, one would actually be 
seeing short pieces of field lines pieced together as they cross the spiral arm and thread through 
the disk. The magnetic field would be mainly azimuthal in between the spiral arms, and only as it 
crosses the arm would it be twisted to align along the arm. Observationally the magnetic field of 
this model would appear the same as that of a large scale magnetic field. 

Similarly, one should see Faraday rotation which is proportional to the amount by which 
the effect of the toroidal field of one sign exceeds the other sign in this region. The amount of 
rotation produced by any one region would be related to difference in areas occupied by the fields 
of different sign. This in turn, is given by the amount by which the initial area at the initial 
position where the radial field is of one sign exceeds the initial area where the radial field is of 
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opposite sign. 

In Parker's argument (Parker 1973a, 1973b) concerning the expulsion of the primordial field 
he introduces the concept of turbulent mixing. Turbulent mixing correctly describes the rate 
at which the mean field will decrease by mixing. However, it can not change the number of 
lines of force threading the disk, since this is fixed by their topology. It only gives the lines a 
displacement. Since, near the edge of the disk, the turbulent motions probably decrease as the 
sources of turbulence do, it is not expected that the lines will be mixed into the halo. Also, 
because of their topology the lines are not lost. Only the length of their extent in the disk can be 
altered by turbulence. Hence, for our model, turbulent diffusion need not destroy the primordial 
field as Parker suggested. Also, blowing bubbles in the lines by cosmic ray pressure will always 
leave the remainder of the line behind, and the total number of lines unaltered. Of course, if the 
lines were entirely horizontal, ambipolar diffusion could destroy the field, since the lines may be 
lifted out of the disk bodily (Parker 1968). 

Therefore, Parker's contention that the primordial field has a short life, need not apply, to 
our model. 

1.1. Our Model 

The model which we consider is very simple. We start with a large ball of plasma whose mass 
is equal to the galactic mass, but whose radius is much larger than the current size of our galaxy 
(figure la). We first assume that the magnetic field filling this sphere is uniform, and makes a 
finite angle a with the rotation axis of the galaxy. 

Then we allow the ball to collapse to a sphere the size of the galaxy (figure lb), and finally 
to a disk the thickness of the galactic disk (figure lc). We assume the first collapse is radial and 
uniform, while the second collapse into the disk is linear and one dimensional. 

During the time when the galaxy contracts uniformly into the disk along the z direction, 
where z is in the direction of f2, we ignore any rotation. Then the resulting magnetic field 
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configuration is as in fig, (lc). The horizontal component of the magnetic field has been amplified 
by the large compressional factor, while the vertical component is unchanged, so that the resulting 
field is nearly parallel to the galactic disk. At this stage, some lines enter the disk from the top 
and leave from the top, e.g. line a. Some lines enter from the bottom and leave through the top, 
e.g. lines b and c. Finally, some lines enter from the bottom and leave through the bottom, e.g. 
line d. It turns out that lines such as a and d are eventually expelled from the disk by ambipolar 
diffusion, so that we ignore them. Now, set the disk into differential rotation at time t = 0. (If we 
include the rotation during collapse the magnetic field will start to wind up earlier. However, the 
result would be the same as though we were to ignore the rotation during collapse, and then after 
the collapse let the disk rotate by an additional amount equal to the amount by which the disk 
rotated during the collapse. This is true provided that the initial field is too weak for ambipolar 
diffusion to be important.) 



1.2. Results of the Model 

We here summarize the conclusions that we found from the detailed analysis of our model, 
given in the body of this paper. 

Initially, after the disk forms, all the lines have a horizontal component that is larger than 
the vertical component by a factor equal to the radius of the disk divided by its thickness 
~ R/ D ~ 100. This is the case if the initial angle a was of order 45 degrees or at least not near 
or 90 degrees. 

Now, because of neglect of rotation in the collapse, the horizontal component of all the field 
lines is in a single direction, the x direction say, so that 

B = Bi± + Bi{Dt&na/R)i = B, t cos Or - B-i sm6§ + B^D tan a/R)z. (1) 

After t = the differential rotation of the disk stretches the radial magnetic field into the toroidal 



direction so that, following a given fluid element whose initial angle is #1, one has 



B = —Bi cos 0i? + Bi(r— t) cos 9i- Bism 9i 9 + Bi(Dtana/R)z. 



(2) 



After a few rotations the second component dominates the first and third components. It is seen 
that the total magnetic field strength grows linearly with time. After the toroidal magnetic field 
becomes strong enough, the magnetic force on the ionized part of the disk forces it through the 
neutral component, primarily in the z direction. 

Consider a single line of force. Let us turn off the differential rotation for a moment. In this 
case, the z motion steepens the line of force, but does not change the vertical field component. 
This leads to a shortening of the line, both in the radial direction and in the azimuthal direction 
(see figure Id). Now, let the differential rotation continue. Because the radial component of the 
magnetic field is reduced, the toroidal component increases more slowly. Eventually, there comes a 
time when the radial field is small enough that the shortening motions in the azimuthal direction 
are stronger than the stretching motion, and the azimuthal component of the magnetic field 
actually decreases even in the presence of differential rotation. From this time on, the magnetic 
field strength decreases at a rate such that the vertical ambipolar velocity is just enough to move 
the plasma a distance approximately equal to the thickness of the galactic disk, in the time t, i. e. 



Now, vo is essentially proportional to the average of the square of the magnetic field strength. 
Therefore, for a uniform partially ionized plasma we have 



where D is the half thickness of the disk, pi is the ion mass density, and where v is the ion neutral 
collision rate. 

Thus, for asymptotically long times one finds that 



v D t = D. 



(3) 



VD ~ 



1 B 2 
PiV 8irD 



(4) 
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That is, the magnetic field strength approaches a saturated time behavior independent of its 
initial value. The saturated time behavior only depends on the ambipolar diffusion properties of 
the interstellar medium, and on the time t. However, the time to reach saturation does depend on 
the initial value of the magnetic field strength. 

The qualitative behavior is shown in figure 2, where the dependence of B on time for different 
initial values is shown. For a very weak initial radial component of the field, saturation is not 
reached during a Hubble time. However, for fields substantially larger than the critical initial field 
strength for reaching saturation, the final saturated field is independent of the initial radial field. 

In the interstellar medium, ambipolar diffusion is not well modeled by diffusion through a 
uniform plasma. In fact, the bulk of the mass of the interstellar medium is in dark clouds, in 
which the degree of ionization is very low. Also, the outward magnetic force is concentrated in the 
volume of the clouds. As a result the ambipolar diffusion velocity is more accurately given by the 
formula 

B 2 (l + f3/a) 

where / is the filling factor for the clouds, pi is the effective ion density in the clouds, v is the 
effective ion- neutral collision rate in the clouds, and [3 /a is the ratio of cosmic ray pressure to 
magnetic pressure. 

The model for the interstellar medium, which we employ to study the magnetic field evolution, 
is sketched in figure 3. The magnetic field is anchored in the clouds whose gravitational mass holds 
the magnetic field in the galactic disk. The magnetic field lines bow up in between the clouds 
pulled outward by magnetic pressure and by cosmic ray pressure. The outward force is balanced 
by ion-neutral friction in the clouds themselves. A derivation of equation 6 based on this model is 
given in section 5. 

Taking plausible values for the properties of the clouds, one finds that the critical initial value 
of the magnetic field, for it to reach saturation in a Hubble time is about 1CT 8 gauss. Further, 
the saturated value of the field at t = 10 10 years is estimated to be 2 microgauss. At this field 
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strength the time to diffuse across the disk D/vj ~ 3 x 10 9 years, which is of order of the lifetime 
of the disk. 

Finally, let us consider the structure of the saturated field. The time evolution of the magnetic 
field refers to the field in the rotating frame. Thus, if we wish the toroidal field at time t and at 
6 = we need to know the initial radial component of the field at 9\ = —Q(r)t. Thus, assuming 
that the field reaches a saturation value of B$, one has 

Be = ±B S , (7) 

where the sign is plus or minus according to the sign of 

[cos(-fi(r)t)]. (8) 

Since Q depends on r, we see that Bg changes sign over a distance Ar such that 
(Afi)t = (Ardn/dr)t « vr, i.e. 

— = ~ n « — (9) 
r r\dn/dr\t \Q\t v ' 

Taking the rotation period of the galaxy to be 200 million years, one finds for t = 10 10 years, 

Ar 100 parsecs. 

Therefore, because the field saturates to a constant field strength, it is predicted that a 
uniform initial field, that has a value greater than > 10~ 8 gauss, immediately after the collapse to 
the disk, leads to a toroidal field that varies with r, at fixed 6, as a square wave with a reversal 
in sign every 100 parsecs or so. Such a field structure would produce no net Faraday rotation and 
would contradict observations. 

However, this result can be traced to our assumption that the initial magnetic field in figure 
1 was entirely uniform. Suppose it were not uniform. Then, after collapse, the radial field is not 
purely sinusoidal (see figure 4). In fact, one expects a behavior more like 

B r (t = 0) = Bi{cos 6 + ecos 26], (10) 

where e is an index of the nonuniformity. As an example, for moderately small e, (cos# + ecos 29) 
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is positive for 



- tt/2 + e < 9 < tt/2 - e, 



(11) 



and negative for 



vr/2-e < 9 < 3vr/2 + e. 



(12) 



Therefore the radial extent in which Bg{t) is negative is larger than the radial extent in which 
it is positive by a factor 



The resulting field is positive over a smaller range of r than the range in which it is negative. 
Consequently, in such a magnetic field there would be a net Faraday rotation of polarized radio 
sources. Note that it is easily possible that the regions in which the radial field is initially weaker 
can end up as regions that dominate in flux over those region that come from regions in which the 
radial field was stronger! Now, if one averages this field over regions much larger than 100 parsecs, 
then the resulting mean field is smooth and axisymmetric. This result is contrary to the generally 
held belief that a primordial field should lead to a field with bisymmetric symmetry(Sofue et al 
1986). This result alone shows that a more careful treatment of the evolution of the primordial 
field, such as that discussed above, leads to quite different results from those commonly assumed. 
(However, as mentioned above, including the effect of the spiral density wave will produce a 
magnetic field lines in the spiral arms that will appear to have bisymmetrical shape.) 

We emphasize that, although the actual magnetic field derived in our model is a tightly 
wound spiral, the magnetic field averaged over a moderate size scale appears to be smooth and 
axisymmetric. 

A second important conclusion of our model is that because of their topology, the field lines 
cannot not be expelled from the disk by ambipolar diffusion. The same line of force diffuses 
downward in the lower part of the disk, and upward in the upper part of the disk. The line must 
thus continue to be threaded through the disk. (See lines b and c in figure lc.) 

A third conclusion that can be drawn from our model is that any single line after saturation 
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has only a finite extent in the disk. For example, if the initial field is a few microgauss, then it 
turns out that the line only extends a radian or so before leaving the disk. This finite extent of 
the magnetic lines would make the escape of cosmic rays from the galaxy possible without the 
necessity of any disconnection of the magnetic field lines in the interstellar medium. 

A final conclusion that should be noted is, that the saturated magnetic field of our model has 
a mean field strength smaller than the rms field strength. This has the consequence that different 
methods for measuring the magnitude of the magnetic field strength should lead to different 
results. The measurement by nonthermal radio emission of the cosmic ray electrons measures the 
rms field strength while Faraday rotation measures the mean strength. 

The outline of the paper is as follows: 

In section 2 an analytic model is developed to demonstrate the properties described in this 
introduction. 

In section 3 a more precise one dimensional numerical simulation is carried out that confirms 
the evolution of the field in the z direction given in section 2. 

In section 4 it is shown that the three dimensional equations for the evolution of the field can 
be reduced to two independent variables. These are z and an angular coordinate u = 6 — Q(r)t 
that is constant along the spirals generated by the differential rotation of the galaxy. A numerical 
simulation of the resulting differential equations is carried out. It is shown that after a long time 
the resulting magnetic field does evolve locally in the essentially same way as is given in sections 2 
and 3. In addition, it varies in radius as a square wave with uneven lobes. 

In section 5 the astrophysics of the interstellar medium clouds is discussed, and a derivation 
is given of expression 6, for the effective mean ambipolar motion of the field in the disk. 

In the concluding section 6, the implications of this model for the evolution of a magnetic 
field of primordial origin are given. The bearing of these implications on the various criticisms of 
the primordial field hypothesis are discussed. 
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2. Local Theory: Analytic 

In this and the next section we wish to consider the local behavior of the magnetic field 
following a fluid element that moves with the galactic rotation. If the ambipolar diffusion were 
strictly in the z direction, the evolution of the field in a given fluid element would be independent of 
its behavior in other fluid elements at different r or 6. It would only be affected by the differential 
rotation of the galaxy, and by ambipolar velocity in the z direction. Thus, in particular, we could 
imagine the magnetic field at different values of 9, behaving in an identical manner. In other 
words, we could replace the general problem by an axisymmetric one. 

Let 

B = B r (r, z)r + B g (r, z)§ + B z (r, z)z. (14) 

Let us neglect the radial velocity, and let the ambipolar velocity be only in the z direction and 
proportional to the z derivative of the magnetic field strength squared. Also, let us neglect any 
turbulent velocity. The only velocities which we consider are the differential rotation of the galaxy, 
Sl(r), and the ambipolar diffusion velocity of the ions. 

Then 

v = rn(r)§ + v z z, (15) 

where 

v. = -K-£-, (16) 

where 

K=(!±fM, (17) 

Pifv 

where as in the introduction (5 /a is the ratio of the cosmic ray pressure to the magnetic pressure, 
Pi is the ion density, and, v is the ion-neutral collision rate in the clouds. The equation for the 
evolution of the magnetic field is 

— = Vx(vxB), (18) 

or, in components, 

()Dr = i-(v z B r ), (19) 



dt dz 
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dB 6 d , . dQ 

ST = + ~dr < 2 °> 

95, <9v 



^2 



dt dr 



B r . (21) 



We integrate these equations numerically in section 3. They are too difficult to handle 
analytically. To treat them approximately, we make the assumption that for all time B r and Bg 
vary parabolically in z i.e. 

B r = B° r (l-z 2 /D 2 ), (22) 
B e = B$(l-z 2 /D 2 ), (23) 



where B® and Bg are functions of time. We apply equation [jl] and equation 2C at z = 0, and 
r = ro, and solve for B®, and Bq as functions of t. 

Thus, making use of equations 16, 18 and 19, we find 

dB r v D 

~dt ~ "d*" (24) 
dB e v D , f dn 

~W = ~~D + {r ^ B - (25) 

v D = ^ B (B 2 r +B 2 ), (26) 

where everything is evaluated at z = 0, r = so we drop the superscripts on B r and Bg. Since 
for galactic rotation rQ, is essentially constant, we have {rdQ,/ 'dr) ro = — O(ro) = — ^o- 

Initially B r and Bg are of the same order of magnitude. Also, we expect Bg to grow, by 
stretching, to a value much greater than its initial value, before ambipolar diffusion becomes 
important. Thus, for simplicity, we make the choice of initial conditions B r = B\,Bg = where 
Bi is the initial value for the radial component of the field. According to equation ||[ B\ = B^ cos 9 
for a fluid element starting at 9, for a fluid element starting at 9. 

Now, if initially Vjj/D <C ^, then according to equation we expect Bg to at first grow 



linearly in time, Then by equation 26, vd increases quadratically in time till Vjj/D ~ S7, and the 



ambipolar velocity starts to affect the evolution of Bg and B r . If the initial field is small, then the 
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contribution of B r to the ambipolar diffusion velocity is never important. In this case we can write 

B 2 

VD = v D1 -£, (27) 



where vdi = KB 2 /(2irD). The solution to the differential equation 24 and equation 25 with vd 
given by equation |27] is 

Br = (1 + ^3)1/2- ( 2§ ) 
Be ~ (1 + ^^3)1/2' ^ 

and 

VD = "D (1 + ^*3) • ( 3 °) 
(The above solution, equations 28 and 29, of equations 24, 25, and 27, is derived explicitly in 
Appendix A. However, it can be shown by direct substitution that it is a solution of equations 24, 
25, and 27.) 



From these equations we see that for 

t <C 



B r is unchanged, Bq increases as £lt, and / vodt <C D, so that up to this time ambipolar diffusion 
carries the plasma only a small fraction of the disk thickness. 



On the other hand, if 

\2v D1 n 2 



( 3D X 1 / 3 



then 



Be ~ Bx i^t = D v¥t = D v (i+l!/a)t' (33) 

and Bq is independent of both the initial value of B r and of fi. It depends only on the ambipolar 
diffusion properties of the interstellar medium. 

For all t, we have 

% - k- < 34 > 
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so for t = 10 years, B r = Bg/300, and the field becomes strongly toroidal. 

The question arises as to how strong B\ must be in order that the saturated solution, 
equation ^3|, is reached. Taking t = tjj, the Hubble time, and making use of the expression for 



vdi, one finds from equation 32 that for saturation to be reached we must have 



Bl > ifi^rm- <35> 

Making use of the numbers derived in section 5 one finds that if Bg is saturated at t = iff, 

then 

D /l0 10 years\ 1/2 /6x 10- 4 \ 1/2 6 , x 

B ° = l^-\ T 1.9 x HT 6 gauss, 36 

lOOpc \ h J \ n i n o J 

where n« is the ion density in the clouds, assumed to be ionized carbon, and no is the mean 
hydrogen density in the interstellar medium. The densities are in cgs units. 

For saturation, the critical value for B\ is 

D /l0 10 years\ 3/2 /6x 10- 4 \ 1/2 „ 

B ler a = — — f 4 x 10- 9 gauss 37 

lOOpc \ t h J \ n i n o J 

where the densities are in cgs units. 

Hence, for the above properties of the clouds and for an initial radial field greater than the 
critical value, the magnetic field at t = tjj saturates at about the presently observed value. Such 
a field arises from compression if, when the galaxy was a sphere of radius 10 kiloparsecs, the 
magnetic field strength was greater than 10 -10 gauss. If before this the virialized radius of the 
protosphere was 100 kiloparsecs, then in this sphere the comoving initial value for the cosmic 
field strength had to be greater than 10~ 12 gauss. (That is if the cosmic field filled all space then 
it had to be so strong that the present value of the magnetic field in intergalactic space must 
now be 10 -12 gauss.) There are good reasons to believe that a magnetic field stronger than this 
minimum value could have been generated by turbulence in the protogalaxy during its collapse to 
this virialized radius (Kulsrud et al. 1996). However, this field would be local to the galaxy and 
not fill all space. 
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Finally, one expects that B z ~ l?i/100, from the initial compression into the disk. Taking B z 
as unaffected by the differential stretching and by the ambipolar diffusion velocity, one can derive 
an expression for the length of a line of force: 



rdd _ B e _ lOOfit / 3D 

~d7-B~ z ~ (i + ^io^3)V2 ~ 100 V 2^? (38) 

or 

. „ / (2.5 x 10" 6 gauss\ , , s 

r A9 = I ^—2 J lOkpc. (39) 

Thus, for example if B\ = 4.5 x 10 gauss, then a line of force passing through the solar 
position in the galaxy, would stretch once around the galaxy. Stronger initial fields lead to shorter 
lines of force. If the initial field is stronger than 2 microgauss, then it is possible that cosmic rays 
can escape along the lines of force into the halo during their average life time. 



3. One Dimensional Numerical Simulation 

It is clear that any gradient of the magnetic field strength will lead to ambipolar diffusion. 
However, gradients in the angular direction are weak, so that we may neglect ambipolar diffusion 
in the angular direction. Similarly, the gradients in the radial direction are weak at first, although 
as pointed out in the introduction, the magnetic field ends up reversing rapidly in the toroidal 
direction, so eventually radial ambipolar diffusion becomes as important as vertical ambipolar 
diffusion. 

In this section, we restrict ourselves to gradients only in the z direction. In this case, the 



axisymmetric approximation is valid, and the relevant equations are equation 1£ and equation 20 



where v z is given by equation 16 



8 B 2 

v z = -K- — . (40) 
oz 8tt 

In the previous section, these equations were reduced to zero dimensions by the ansatz that 



B r and B z were parabolic in z according to equation \22 and equation 23, and the basic equations 
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were applied only at z = 0. In the present section we treat these equations numerically for all z 
with \z\ < D, and drop the parabolic assumption. We treat the disk as uniform so that pj and v 
are taken as constants. 

We need boundary conditions at z = ±D. We expect B to be quite small outside the disk, 
\z\ > D. (This is because we suppose that neutrals are absent in the halo and u = 0, so that the 
flow velocity v z becomes very large. Since flux is conserved, v B must be a constant in a steady 
state and the magnetic field must be very small.) In order to match smoothly to the outer region, 
we first assume that pi is a constant for all z and is very small. In addition, we assume that v is 
constant for \z\ < D, and decreases rapidly to zero in a narrow region, D < \z\ < D + A. Then 
because v in the halo region, \z\ > D, v z becomes so large that inertia is important. 

The equation for v z should read 

(dv z dv z \ dB 2 /8n , 

where we have set piv' = (piU / f)(\ + j3/a)/f. Because the time evolution is over the Hubble time 
tn, dv z jdt <C v z dv z /dz, i.e. v z /D is very large compared to 1/tu- Thus, we drop the partial time 
derivative. 

The equations for B r and Bg are the same as in the disk. However, because v z is large the 
divergence terms J^(v z B r ) and J^(v z Bg) are much larger than the other terms in the region near 
\z\ = D, and in the halo. Thus, we have 

Yz (v z B)=0, (42) 



where B = y B% + Bg is the magnitude of the magnetic field. This means that 

vB = $ = const, (43) 

where $ is the rate of flow of magnetic flux, and it is a constant in these regions. (We drop the 
subscript on v) 

Dividing equation 41 by pi and dropping the partial time derivative, we get 

l(" 2 + 4)^ 2,/ - <44> 
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Combining this equation with equation ^ we have 

d_( 2 $ 2 ' 

dz \ A^piV 12 



-2i/v, 



or 



or 



(- 



$2 



1 



v'dz. 



(45) 



(46) 



(47) 



l 47T/9if 4 

For z larger than D + A, the right hand side becomes a constant, so that v A — > <I )2 /47rpj. As 
z decreases below the right hand side becomes linear in z, and we find that for z far enough 
into the disk, the inertial term becomes negligible. Hence, 

$2 



3V 3 



(48) 



and v approaches a small value. This result breaks down when v is small enough that the other 



terms in equation 19 and equation 20 become important. 



Thus, the connection of the main part of the disk to the halo is given by equation 47 

/ $2 



1 dv 



v' 'dz, 



(49) 



where v c = ($ 2 /4-7rpi) 1//4 . Transforming this to an equation for B we have 



B / B 4 



$>dB 



v'dz, 



(50) 



l Bc \±1Tpi$ 2 ) B 2 

where B c = &/v c = (47r/jj<I> 2 )~ 1 ' /4 , and where z c is the value of z where B = B c . Let us write 
v'dz = v' A', where on the right hand side v' denotes the constant value of v' in the disk. Then 



for z far enough into the disk, equation [5C] reduces to 

B 3 



(D - e - z)v Q , 



where 



e = A' + 2v c /3u Q . 



(51) 



(52) 
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Now, A' is small by assumption. If we estimate <3? as -Bo^effj where v oeti 
is the order of the ambipolar diffusion velocity at the center of the the disk, then 
v c ~ (Bov Dg q jyj Airpi) 1 / 2 <C Bo/y/Airpi, so that e is much smaller than the distance an 
alfven wave can propagate in the ion-neutral collision time. This is clearly a microscopic distance 
compared to the thickness of the disk, so that we may neglect it. 



In summary, equation 50 implies that the inner solution essentially vanishes at z ~ iD, so 
that we may take our boundary condition to be B r = Bq = at z = D. Given this solution the 
above analysis shows how to smoothly continue it into the halo. 



Equation 1£ and equation 20 can be made dimensionless by a proper transformation. We 
choose a transformation that is consistent with that employed in the next section. Define a unit 
of time to by 

fVo = R/D, (53) 

where R is the radius of the sun's galactic orbit, and f2o = 0(i?) is its angular velocity. Next, 
choose a unit of magnetic field Bq to satisfy 

M = £. (54) 
AttD t 

For such a field the ambipolar diffusion velocity is such that the field lines cross the disk in a time 
to- Now, let 

t = t t', (55) 
B r = (D/R)B B' r , 

Bg = B B' e , 

z = Dz'. 

For a cloud ion density of 6 x 10 _4 /cm 3 , and a mean interstellar medium density of 
no = 10/cm 3 ,5o = 2.85 x 10 _6 gauss, (see section 5). Also, D = 100 parsecs, and to = 3 x 10 9 
years. The Hubble time in these dimensionless units is about 3. 
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The details of the numerical simulation are given in another paper, (Howard, 1996). The 
dimensionless equations to be solved are 

dB' r 1 d ( dB'\ 



dB' a 1 d ,dB 



12 



dt> 2dz' y d 



l^ B 'e)- B r- (57) 



Our dimensionless boundary condition is B' r = BL = at z = ±1, 

The results of the integration of these equations are shown in figures 5 and 6. The initial 
profiles of B r and Bq are parabolic. 

In figure 5 the initial value of B' T is B; L = 0.01, so that B r = 2.85 x 10~ 8 gauss. It is seen that 
the qualitative behavior is the same as that described in sections 1 and 2. The field at first grows 
and then relaxes back to a decaying solution. Figure 6 gives the time evolution of B' e at z = 
for the set of initial conditions, B r = Bi cos theta where 6 runs from zero to one hundred eight 
degrees in increments of fifteen degrees. It is seen that the curves all tend to the same asymptotic 
behavior and are similar to those in figure 2. As in figure 2, the field from weaker initial values of 
B r takes longer times to reach its peak value and the saturation curve. These results are similar 
to those which are represented by equations 28 and 29. 

It is seen in figure 5 that, except near the edge, the profile remains similar to a 
parabola. Near the edge the cube root behavior of equation ^l] is also evident. This cube 
root behavior can be understood since the flux v z B is roughly constant, and therefore 
(dB' 2 /dz')B' » (0/0*0(1 - z'fl 3 x (1 - z') 1 / 3 w (1 - z')° = a constant. 

The most important feature in figure 5 and 6, is the saturation of B to the envelope curve. 
This saturation occurs if the initial value of B' r is larger than 0.01. Thus, even in our more precise 
calculations, information on the initial B is lost, and the final value of B at fixed time depends 
only on the initial sign of B r . 
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4. Two Dimensional Analysis 

In the last section we have treated the evolution of the magnetic field under the influence 
of differential rotation and ambipolar diffusion in a one dimensional approximation. Only the z 
component of the ambipolar diffusion was kept. In this approximation the magnetic field in each 
column of fluid at r, 9 evolves independently of any other r, 6 column of fluid. 

In early times, because the galactic disk is thin compared to its radius, this is a good 
approximation, since the horizontal gradients of B 2 are small compared to the vertical gradients. 
However, because of the differential rotation of the galaxy, fluid elements at different radii, that 
were initially very far from each other are brought much closer together, and the horizontal 
gradients are increased until they become comparable to the vertical gradients. 

For example, two fluid elements that were initially on opposite sides of the galaxy, and at a 
difference in radius of one hundred parsecs, will be brought to a position on the same radius after 
about fifty rotations. Since the evolution of the field in these two fluid elements is very different, 
we expect the gradient of B 2 in the radial direction to become as large as that in the vertical 
direction. As a consequence, ambipolar drift velocities in the horizontal direction can be expected 
to be as large as those in the vertical direction. On the other hand, two fluid elements at the same 
radius, but initially far apart, remain far apart. The ambipolar motions in the 6 direction should 
remain small. 

Thus, we expect that, at first, only ambipolar z motions are important for the evolution of 
the magnetic field, but that eventually the radial ambipolar motions also will become important, 
although not the angular ambipolar motions. That is, we expect the problem to become two 
dimensional. 

In order to properly demonstrate this evolution, we introduce a new independent variable 

u = 6- Q(r)t. (58) 

If ambipolar diffusion is neglected, this variable is just the initial angular position of a fluid 
element that is at the position r, 6, z at time t. Since in the absence of ambipolar diffusion the 
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variables r, it, and z are constant, following a given fluid element. They would be the Lagrangian 
variables if only rotational motion is considered. It is appropriate to describe the evolution of the 
magnetic field components B r ,Bg, and B z in terms of these variables. 



Inspection of equation 58 shows that when Qt is large, u varies rapidly with r at fixed 0, in 
agreement with the above qualitative discussion. Changing r by only a small amount will change 
the initial angular position by n. Thus, we expect that the components of B will vary finitely 
with u, but only slowly with r for fixed u. The surfaces of constant u are tightly wrapped spirals. 
Thus, the behavior of the field should be finite in r, u, and z, but only gradients with respect to u 
and z should be important. 

To see this, let us first write the total velocity, v = w + Qr#, where w is the ambipolar 
velocity. Next, let us derive the equations for the components B r ,Bg,B z and w r ,wg,w z in terms 
of the Eulerian variables r, 9 and z. 

^ = (B • V)w r - (w • V)B r - B r (V • w) - fi^, (59) 
dB e B e w r B r w e dB r 

-—- = (B • V)W0 (w • VjBe H ^^Z" ~ B e(y ' w )> ( 60 ) 

at r r oo 

^ = (B • V)w z - (w • V)B Z - - B Z (V • w), (61) 

where the ambipolar velocities are: 

8B 2 /8Tr 8B 2 /%tt T dB 2 /8n . . 

w r = -K J ,wg = -K ^-,w z = -K J , 62 

or rot) oz 



and 



V w = — + — + — + — (63) 
dr rdO r dz 



Finally, let us transform these equations to the new coordinates r, u, z. In doing this we assume 
that the galactic rotation velocity v c = VLr is a constant. 

The result of the transformation is 

dB r = B^dw L _^ B dwr , g4 , 

dt r du 2 dz 

dB r wg dB r w r dB r dB r 

— W r ^ ^ ill r ID, 



dr r du r du dz 
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r , dwe_ dw z 

'' r rdu dz 

9B e = B dwe_ SU B dwe_ 

dt r dr r r du 

„ dwg dBg Wq dBg fit dBg 

+B z~^ w r~Z a W r-Z— 

dz or r du r du 

dBe 

Oz 



-w z — QB 



dw r Qt dw r dw z 
-Bq-z B e — Bg——, 

or r ou oz 

9B Z _ B dw z | B e dw z | ^B r dw z 

dt r dr r du r du 

dB z w e dB z w r dB z 

~ w r~^ 3 nt 5 — 

dr r du r du 

dB z dw r 

~ W z-5> B z~^~ 

dz dr 

B z dw r D w r B z dw e 

-lit t> z — , 

r du r r du 

dB 2 /8ir ntdB 2 /8ir 

w r = -K—£ K — , (67) 

dr r du 

KdB 2 /8iT , , 

w e = 68 

r du 

w z = ~K—£-. (69) 

Only a few of these terms are important. To see this let us introduce dimensionless variables 
for the velocity and field components. (We will denote the dimensionless variables by primes.) We 
will choose these variables based on our one-dimensional results, as follows: 

The unit of length for the z variable is the galactic disk thickness D. 

z = Dz'. (70) 

The variation of quantities with r is finite over the distance R the radius of the sun's orbit in the 
galaxy, so we set 

r = Rr'. (71) 
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The variable u is already dimensionless and quantities vary finitely with it, so we leave it 
unchanged. 

The unit of time to should be of order of the age of the disk. During this time the number of 
radians through which the galaxy rotates, fit is of the same order of magnitude as the ratio R/D, 
so for analytic convenience we choose to so that 

fVo = R/D, (72) 

and set 

t = t t'. (73) 
If we take R/D = 100, and Q = 2tt/(2 x 10 8 years), then t = 3 x 10 9 years. 

It is natural to choose the unit for Bg, Bq as that field whose z gradient produces an ambipolar 
z velocity of order D/to, that is an average velocity near that which would be produced by the 
saturated field. Thus, we choose Bo so that 

KB$/(4irD) = D/t , (74) 

and set 

Be = B B' e . (75) 
(Note that the definitions in these units are consistent with those of section 3.) 

In most cases of interest to us, B r and B z are much smaller than Bo. (In general Bg is 
initially small compared to Bo, but it grows by a factor of R/D to up to the saturated value of 
about Bo-) Thus, we set 

B r = (D/R)B B' r . (76) 

The vertical field B z starts out even weaker than this magnitude since the initial horizontal 
components of the magnetic field were amplified by the initial compression which formed the disk. 
However, the B z field is amplified up to the size of the B r field by the shear of the radial ambipolar 
velocity acting on B r . We thus transform B z by 

B z = (D/R)B B' Z . (77) 
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The ambipolar velocities w r ,wg, and w z arise from gradients of the magnetic field in 
the corresponding directions. Thus, after Bg has been amplified by stretching to be of order 
Bq,w z ~ KBq/D. Similarly, after the differential rotation has acted to reduce the scale of the 
variation of B 2 in the radial direction, w r becomes of order w z . However, the scale of variation 
in the 9 direction remains of order R over the age of the galaxy so that wg « KB 2 /R. Thus, we 
change the w components to w' components by 

w r = (D/toK, (78) 
w e = (D 2 /Rt )w' e , 
w z = (D/t )w' z . 

Now if we transform equations 64 to 69 by the change of variables equations 70 to 78, and 
clear the dimensional factors to, Bo etc. from the left hand side, we find that the terms on the 
right hand side are either independent of dimensional units entirely, or are proportional to powers 
of (D/R) « 1. The full equations are given in appendix B. Dropping these "smaller" terms 
proportional to a power of D/R greater than zero, we find that the equations for the dimensionless 
variables, to lowest order in D/R are 

M 1 ,, f dB> r d(w' z B' r ) 

3t> r' 2 r 8u dz> { ' 

, B 'e dw 'r , B l dw 'r 

r' du z dz' ' 
dB' e t' d{w' r B' e ) dw' z B' g 



dt' r' 2 du dz' 

B' r 



(80) 



dB' B'ndw' t'B'dw' 

of r ou r z ou 

t'w' r dB' z ,dB' z t' .dw'r 

W~— -^r B . 



r' 2 du z dz r' 2 z du 

t' dB' 2 , , 

W '=-2^^ (82) 
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< = -2~37- (84) 

Note that w' e does not occur in the equations for the evolution of the B' components. Also, the 
r' derivatives are absent from these lower order equations. The initial conditions on B'(t',u, z';r') 
are 

B' r (0,u,z',r') = B r (0,r'R,u,z'D), (85) 
B e (0,u,z',r') = B e (0,r'R,u,z'D), 
B' z (0,u,z',r') = B z (0,r'R,u,z'D). 

These transformations are formal, but they enable us to correctly drop the terms whose effect 
is small. Once these terms are dropped, the equations reduce to two dimensional equations, which 
are more easily handled numerically. Although we have assumed that the dimensionless variables 
are originally of order unity, it may be the case that they differ substantially from unity. However, 
an examination of the various possible relevant cases leads to the conviction that all the important 
terms have been kept as well as other terms which are, perhaps, unimportant. For example, the 
initial value of B' d is much smaller than unity. However, because of the shearing terms in equation 
80 (the last term) B' e grows to finite order when t' becomes of order unity, so during the later 
stages of the galactic disk, B' e is of order unity. 

Many of the terms in equations 79 to 81 have an obvious significance. The second term 
on the right hand side of the B' equation is the vertical decompression term present in the one 
dimensional simulation. The first term is the radial decompression term, which only becomes 
important when t' « 1, and the wrapping up has made the radial ambipolar diffusion important. 

Similarly, there is a z decompressional term in the B' r equation, but, of course, no radial 
decompression term. There is a term representing the effect of shear on the toroidal field in 
increasing the radial component, and a similar term resulting from the action of shear on the 
B' z component. These shear terms would be small if Bg were of order of B r , or if B z were much 
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smaller than B r , which is the case initially. However, B z is increased by shear terms over the age 
of the disk, to a value considerably larger than its initial value. 

Equations 79 to 84 only contain derivatives with respect to t', z' and u, and none with respect 
to r'. Thus, r' is only a parameter in these differential equations. Therefore, the components of 
the magnetic field evolve independently of those at a different value of r' (to lowest order in D/R) . 

Further, u does not occur explicitly in the differential equations. The initial conditions, 
equations 85 do involve u, and are periodic in it. Therefore, the magnetic field components remain 
periodic in u for all t'. 

Let us consider the behavior of such a solution, periodic in u, in the neighborhood of the sun, 
r' = 1, at fixed z', say z' = 0. Transform the solution back to r, 9 coordinates. For fixed r and t, 
the solution is periodic in 9, e.g. Bg(r,9) = Bg(r,u + Q,(r)t). Moreover, for fixed 9 we can write 

Ar 

u = 9-Q(r)t = 9 Q t (86) 

= 9-^ = 9-^ 
r D D 

so 

Ar = 9-uy. (87) 

Thus, for fixed 9, r changes by an amount 2irD/t' when u changes by its periodic length 2ir. 
Since r' changes by a small amount f» 2ir D/Rt', we may ignore the dependence of the solution 
for the components of B on r' and the components of B are nearly periodic in r ( at fixed 9 ). 
However, because of the actual dependence of the solution on r', as a parameter in the equations, 
the amplitude and phase ( as well as the shape) of the periodic solution do change slowly when 
one goes a distance comparable with the radius of the galaxy. 

Equations 79 to 84 were integrated numerically. The details of the integration are discussed 
in Howard(f995, f996), where most of the results are presented. Initial conditions were set by 
starting with a cosmic field before compression into the disk, and then calculating the resultant 
fields. In this paper we present the results for two initial cases. Only the results for the integrations 
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at r = R, the radius of the galactic solar orbit, are included. The variation of Bg as a function 
of u, at z = 0, is plotted in figure 7 for the case that the initial cosmic field was uniform. The 
antisymmetry in u is evident, and it is clear, after transforming the field to be a function of r as 
the independent variable by equation 87, that no Faraday rotation would be produced by this 
field. The same result for the case when the initial cosmic field was nonuniform is shown in figure 
8. [The initial cosmic field was chosen so after compression into the disk the horizontal field was 
B = Bi[.5 + x)x + yy).] The resulting saturated field is not antisymmetric, and does not average 
out in u. It also would not average out when transformed to be a function of r, and would produce 
a Faraday rotation. The variation of Bg with z at u = at a time of 9 gigayears, is shown in 
figure 9. It has the parabolic shape found in section 3. The variation of Bg with time at the point 
u = 0, z = 0, is given in figure 10 for the two cases. The results are also similar to those of figure 2 
derived from the simple parabolic approximation. 

5. Ambipolar Diffusion in the Interstellar Medium 

We now consider the averaged equations for the magnetic field, taking into account the 
interstellar clouds. The bulk of interstellar matter is in the form of diffuse clouds and molecular 
clouds. Because the properties of the molecular clouds are not very well known, we make the 
simplifying assumption that essentially all the interstellar matter is in diffuse clouds, with a small 
amount of matter in the intercloud region. In describing the clouds, we make use of properties 
given by Spitzer(1968). 

We assume that all the clouds are identical. We further include the cosmic ray pressure, and 
the magnetic pressure, in the intercloud region, but neglect the pressure of the intercloud matter. 
Then the cosmic rays and the magnetic fields are held in the disk against their outward pressures 
by the weight of the clouds in the gravitational field of the stars. (See figure 3.) 

Now the force due to the magnetic and cosmic ray pressure gradients is exerted only on the 
ionized matter in the clouds, while the gravitational force is exerted mainly on the neutrals in the 
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clouds, since the fraction of ionization in the clouds is generally very low. Thus, these contrary 
forces pull the ions through the neutrals with ambipolar diffusion velocity vp. The frictional force 
between the ions and neutrals is proportional to vp. By equating the magnetic plus the cosmic 
ray force to the frictional force, we can obtain the mean ambipolar velocity in the clouds. Now, we 
assume that the cosmic ray pressure pr is related to the magnetic pressure B 2 /8ir by the factor 
/3/a(Spitzer 1968). We take (3/a independent of time and space. This is plausible since when 
the magnetic field is strong, we expect the cosmic ray confinement to be better and therefore the 
cosmic ray pressure to be larger. 

The mean vertical force per unit volume produced by the magnetic field strength gradients 
and cosmic ray pressure gradients is 

+ ,/„)**£!. (88) 

But this force is counterbalanced by the gravitational force on neutrals in the clouds, which occupy 
a fractional volume equal to the filling factor, /, times the total volume. 

Thus, the force per unit volume on the ions in the clouds is 

F doud = j = -(l + p/ a )j^^. (89) 

This force produces an ambipolar velocity, vd, of the ions relative to the neutrals such that 

Fdoud = riim*uv D , (90) 

where m* is the mean ion mass, v is the effective ion-neutral collision rate for momentum transfer 
and Hi is the ion number density, in the clouds. Now, 

ttih 

v pa n c < av > (91) 

m* + rriH 

where a is the momentum transfer collision cross section, and we assume that the neutrals are all 
hydrogen with atomic mass m#. Thus, we have 

777# 777* 

Fdoud = rii— — < av > n c v D (92) 

m* + run 
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If m* >> run, then the mass factor is ra#. Hence, if the ions are mostly singly charged 
carbon, the mass factor is « run, while if they are mostly hydrogen ions then it is m#/2. The 
ions are tied to the magnetic field lines, since the plasma is effectively infinitely conducting. If 
there are several species of ions, they all have the same cross field ambipolar velocity. 

Thus, the ambipolar diffusion velocity in the clouds is 

l (l + /?/q)V(i? 2 /87r) 

v d = 7 ; , (93) 

/ riiinlff < av > n c 

where Ttg// =< m*mH /{rn* + tuh) > is the effective mass averaged over ion species, and a is the 
ion-neutral cross section. 

The velocity vd is the mean velocity of an ion inside a cloud. It is also the velocity of a 
given line of force. Now, the ions are continually recombining and being replaced by other ions, so 
it is actually the motion of the magnetic lines of force that has significance. Further, the clouds 
themselves have a short life compared to the age of the disk. They collide with other clouds every 
10 years or so, and then quickly reform. After the collision the cloud material is dispersed, but 
because of its high conductivity, it stays connected to the same lines of force. After the cloud 
reforms the lines of force are still connected to the same mass. The lines then continue to move in 
the reformed cloud, again in the opposite direction to the field gradient. As a consequence, during 
each cloud lifetime, the lines of force pass through a certain amount of mass before the cloud is 
destroyed. 

The amount of mass passed during a cloud life time, by a single flux tube of diameter d is 
2v£>t c adp c where t c is the cloud life time, 2a is the cloud diameter, and p c is the cloud density. 
However, because of our assumption that the bulk of the matter is in these diffuse clouds, it 
must be the case that there is only a short time in between the destruction of one cloud, and 
its reformation. Thus, during a time t, the field passes through « t/t c successive clouds. If we 
consider a length of a line of force L, at any one time it passes through fL/2a different clouds. 
Then, in the time t, the amount of matter passed by this length of a given tube of force (of 
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diameter d ) is 



AM = 2adv D t cPc —^ = Lfp c dv D t. (94) 
la t r 



But fp c , is the mean density, p, of interstellar matter. Thus, if we define Ax as the average 
effective distance that this magnetic tube moves through the disk by 

AM = pAxLd, (95) 

we get 

dpAxL = Ldpvot, (96) 

or 

^ = vd- (97) 

This effective velocity, averaged over many cloud lifetimes, is the only velocity for the magnetic 
field lines that makes sense. It is equal to the velocity, vd, of the ions in the cloud material. 

So far, we have made the simplifying assumption that the clouds are stationary, which is of 
course not the case. As the clouds move through the interstellar medium they stretch the lines, 
and additional tension forces and ambipolar velocities arise. However, these velocities are always 
directed toward the mean position of the cloud material, and thus they tend to average to zero. 
The ambipolar velocity which we have calculated above, under the assumption of stationary 
clouds, actually gives the rate of displacement of the lines of force relative to the mean position of 
the cloud. It is a secular velocity, and it is the only velocity which really counts. 

We assume that in the diffuse interstellar clouds only carbon is ionized, so that m* e ff = run- 
We choose < av >= 2 x 10 _9 cm 2 /sec(Spitzer 1978). If the filling factor of the clouds / = 0.1, 
and if the mean interstellar density of hydrogen is no ~ 1.0/cm 3 , then the density in the clouds 
is n c ~ 10/cm 3 . The cosmic abundance of carbon is 3 x 10~ 4 (Allen 1963). We assume that this 
abundance is depleted by Co ~ 0.2, (Spitzer 1978), the rest of the carbon being locked up in grains. 
We further take j3/a = 2. Then comparing equation ^ with equation ^ gives 

K = (1 + /3/a) = 1.5 x 10 36 , (98) 

riim* e ff < av > uq 



in cgs units. Taking D = 100 pc, we get from equation 54 

B = 2.85 x 10~ 6 gauss. 
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We take t = 3 x 10 9 years. We see that if t' = 3, corresponding to an age for the galactic disk 
of 9 x 10 9 years, then the present value of the magnetic field from figure 8, should be 1.5 x 10 -6 
gauss. 

This value depends on our assumptions concerning the properties of the clouds. In particular, 
if the hydrogen in the interstellar clouds is partly ionized by low energy cosmic rays penetrating 
the clouds, then the ambipolar diffusion will be slower, and K will be smaller, leading to a larger 
value for Bq. The consequence of this is, that the initially dimensionless magnetic field B' r will be 
smaller, so that it would take longer to reach saturation. However, the saturated value will be 
larger. 

6. Conclusion 

We have assumed that there was a cosmic magnetic field present before the galaxy formed. 
On the basis of this hypothesis we have constructed a simplified model of the galactic disk in order 
to investigate how the magnetic field would evolve, and what it would look like at present. The 
two essential ingredients of this model are: the differential rotation of the interstellar medium, and 
the motion of the field lines through the interstellar medium produced by the ambipolar diffusion 
of the ionized component of the plasma driven through the neutrals by magnetic pressure and 
cosmic ray pressure gradients. The effect of turbulent motions is assumed to average out. No large 
scale mean field dynamo action on the magnetic field was included in this model. (This model 
differs from that of Piddington in that the magnetic field is too weak to effect the galactic rotation 
of the interstellar medium, and also, ambipolar diffusion is included.) 

The consequences of the model were investigated by an approximate analytical calculation 
in Section 2, and by more detailed numerical simulations in Section 3 and 4. These simulations 
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confirmed the results of the approximate analysis of section 2. 
The basic results are: 

(1) To first approximation, the magnetic field evolves locally following a rotating fluid 
element. It first grows by stretching the radial component of the magnetic field into the toroidal 
direction. When the field becomes strong enough, the line commences to shorten because of 
the vertical motions produced by ambipolar diffusion. This reduces the radial component, and 
therefore the stretching. After a certain time, the field strength saturates and starts to decrease as 
the reciprocal square root of time. This asymptotic behavior is determined only by the ambipolar 
diffusion properties of the clouds. Thus, the field strength everywhere approaches the same value 
at a given time. At the present time, this value is estimated to be in the range of a few microgauss. 
This saturated value is independent of the initial value which the field had when the disk first 
formed, provided that the initial value of the field strength is greater than 1CT 8 gauss. The extent 
of each magnetic field line in the toroidal direction also saturates, but its length in the disk does 
depend on the initial value. 

(2) The direction of the toroidal field in any given fluid element depends on the sign of the 
initial radial component. Since this sign varies with position, and since differential rotation mixes 
these positions, it turns out that the resulting toroidal field varies rapidly with radius along a fixed 
radial direction. The toroidal field changes direction on a scale of a hundred parsecs. Because the 
saturated field strength is nearly constant in magnitude, the toroidal field strength as a function 
of r at fixed 6 varies as a square wave. However, the lobes of this square wave need not be equal 
since the regions of one initial sign of B r may be larger than those of the other sign. In this case, 
the model predicts a toroidal field that would produce a net Faraday rotation in radio sources 
such as pulsars or polarized extragalactic radio sources, in spite of its rapid variation in sign. 

It is the prevailing belief that the galactic magnetic field does not reverse in radius on 
small scales. In fact, this belief is grounded in an analysis of Faraday rotation measures of 
pulsars (Hamilton and Lynn 1987). In analyzing these rotation measure Rand and Kulkarni(1988) 
employed various simple models of the galactic field. In every one of these models, the magnetic 
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field reversed only on large scales of order one kiloparsec. These models, which only allowed a slow 
variation of B, led to results that were consistent with the observations, and thus supported the 
general belief in reversal only on large scales. However, when this analysis was carried out, there 
was no apparent reason reason not to consider a model in which the field varied rapidly in radius 
on scales of order of a hundred parsecs, although in hindsight it could have been done. However, 
such a model was not included in the analysis of the rotation measures, and thus it was not tested. 
Therefore, at present, there is no reason to exclude such a model. 

In short, the prevailing belief that the galactic magnetic field is of constant sign on a large 
scale, actually resulted from the assumption that the field was a large scale field, and from the 
consistency of this assumed model with observations. Therefore, this conclusion has not been 
rigorously demonstrated. A magnetic field, such as that arrived at from our model, which has the 
rapid variation of the field with radius, would lead to fluctuations from the mean. Indeed, such 
fluctuations are in the data and are attributed to a general isotropic random magnetic field in 
addition to the mean field. 

We feel the the non uniform square wave could probably fit the observations equally as well 
as the other models, so that our model can also be shown to be consistent with observations. We 
have not yet demonstrated quantitative consistency with observations, but we hope to carry out 
this task in the future. 

(3) The magnetic field observed in our galaxy and in other galaxies should actually be the 
average of the true detailed magnetic field averaged over regions in space larger than those regions 
over which we find our field to vary (several hundred parsecs at least). Thus, under averaging 
the field of our model would actually appear as an axisymmetric toroidal magnetic field. An 
analysis of the various galactic magnetic fields has made the ansatz that the origin of the field 
can be distinguished as primordial if it has bisymmetric symmetry, and as due to a dynamo if it 
has axisymmetry(Sofue et al. 1990). (It must be borne in mind that the actual magnetic field is 
perturbed by the spiral arms, so that it is parallel to the arms in the region of the spiral arms and 
toroidal in the region in between the arms as described in the introduction.) On the basis of our 
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analysis, we conclude that this method of distinguishing the origin is not valid, and could lead to 
incorrect conclusions concerning the origin of the galactic magnetic fields. 

4. Parker(1968, 1973a,b) and others (Ruzmaikin et al 1988) have put forth three objections to 
a primordial origin. These objections are: (a) A primordial field would be tightly wrapped up, 
contrary to observations, (b) Such a field would be expelled by ambipolar diffusion or turbulent 
diffusion, (c) There is no known mechanism to produce a large scale primordial field in the early 
universe. 

On the basis of our model we can counter the first two objections (a) and (b). Our model 
does indeed lead to a tightly wrapped spiral magnetic field. However, when averaged over a 
sufficiently large scale, as is done automatically by observations, the resulting field should actually 
appear to be axisymmetric and azimuthal. This averaged field is, thus, not in disagreement with 
observations. In addition, if the large scale field is not uniform, then the field after saturation by 
ambipolar diffusion creates larger regions of one sign then of the other sign and the averaged field 
is not zero. Thus, objection (a) does not defeat our model for a field of primordial origin. 

With respect to objection (b) our model predicts field lines which thread through the disk, 
entering on one edge and leaving on the other. Thus, it is impossible for a vertical ambipolar 
motion or turbulent mixing to expel lines of force. This result counters objection (b). 

The third objection is still open to debate and we do not discuss it. 

The two observations that should test our model are: The magnetic field in other galaxies 
should be observed to be toroidal in between the arms. A reanalysis of the pulsar rotation 
measures should fit our model without very large extra fluctuations. These fluctuations should be 
accounted for by the rapid reversals of the toroidal field predicted by our model. 

To summarize, we have shown that a careful analysis of the evolution of a primordial field 
throws new light on the way one should view the primordial field hypothesis. 
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A. Appendix A 

In this appendix we give the derivation of the solution given in equations 28-30, to the 
differential equations (24), (25), and (27) of section 2. Let us drop the 6 subscript in Bg. Then 
the equations to be solved are 

^ = - V -^B r (Al) 
dt D y ' 

_ = - v -R B -n Br 

V 4v D = kB\ 



D 



Define £ by- 



Then 



or 



and 



d T = ^ = kB 2 . (A2) 
dt D v ; 



f = -B r , (A3) 
B r = B ie ^, (A4) 



!{ = - b -up = - b -W^ (A5) 



where C = VLB\/k is a constant. 
Multiplying by e^, we get 

dC~^ ' ~ B 2 ~ {Bet) 2 
Now, B is negative. Let y = —Be^. Then 

dy _ Ce 2 Z 

Integrating and assuming the initial value of y is small we get 



^(BeZ) = — = —j-,. (A6) 



(A7) 



^ = f(e 2? -l), (A8) 
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or 



,3C 



(^i)l/3( e 2f _ 1)V3. 



(A9) 



Now, 



e 2^ = kB 2 e 2^ = k 2 

dr y 



(A10) 



Next, let = rj. Then 



^ = f(f) 2/3 (,-D 2/3 , 



(All) 



or 



Integrating this equation we have 



(A12) 



(A13) 



Now expanding this equation and making use of the definition of r] we have 



e 2 « = 1 + -C 2 k 3 t 3 t 3 . 
3 



(A14) 



The toroidal field B = ye * , and y is given by equation |A9| so 



D 2 2-2* ,3C. 2/3 (e 2 «-l) 2 / 3 
B 2 = y 2 e ^ = (—Y'^ ^ — ' 



(A15) 



or 



C 2 kH 2 



1 + (2/3)C 2 £; 3 t 3 ' 

Restoring the definition C, and setting vdi = fe-B? we get equation (29) of section II for Bq 
Equation (28) for B T is obtained from equation A4 and equation A14| . 



(A16) 
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B. Appendix B 

In this appendix we give the full equations for dB' r /dt,dB' g /dt and dB' z /dt including all 
the terms in R/D. They result from the transformation of equations 64 to 69 to dimensionless 
variables by making use of equations 70 to 77. 

For convenience, the terms are listed in exactly the same order as in the original dimensional 
equations. The resulting equations are 

OK B^dw^_ M 

dt' r> dz' z dz' { ' 

£>\ , dB' r m\ 2 w' 9 dB' r w' r dB' r dB' r 
R J r dr' \RJ r' du r' 2 du z dz' 
(£\ < B '- (D\ 2 Kdw^_ B ,dw' z 



\RJ r' r \R r' du r dz 



dB e 
dt' 



+ 



—w 



RJ 

RJ 

dz' 
D 

~R 



r Q r , 

t dw' e 
dz' 
BL 



D\ 2 B^dw^ t , 



RJ 



- — w: 



,,/2 



du 



j dB e 
dr' 



B, 



i dw' r t' , dw' r , dw' z 



dr' 



dr' 



-B' a 



dz' 



w, 



dB' e t' , dB' e dB e 



du 



-w 

r' 2 r du dr' 



(B2) 



dBi 

dt' 



£\ ,dK 

R r dr' 



+ 



B' a dw'. 



t' dw' 



fD\ ,dB' z 

~ (r) w ^ 



du 

D\ 2 w' g dB' z 



du 



(B3) 



-w. 



,dBi 

dz' 



r' 2 2 du 



R 



RJ r' 
, dw' r 

~dV 

B' z , 
—v-wL - 



du 



t' , dB' z 
— w - 

r' 2 r du 



2 B' dw'n 



du 



The equations for the components of w are: 



w' g 



t' dB' 2 dB' 

w = 1 

r' 2 du dr' 

i dB' 2 



dB 



12 



dz> 



In these equations B' 2 stands for 



B'" = B' 2 + (^]\b 2 + B 2 ) 



RJ 

We do not expand the expression for the components of w' out fully. 
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FIGURE CAPTIONS 

Fig. 1 (a) The protogalaxy before initial collapse, showing a uniform cosmic field threading it. 
The magnetic field makes an angle a with the rotation axis f2. (b) The magnetic field after 
the protogalaxy has collapsed to the present size of the galaxy but before the galactic disk 
has formed. The lines remain threaded through the galactic plasma and connected to the 
rest of the cosmos, (c) The magnetic field lines after collapse to the galactic disk. Some 
of the lines such as a and d will escape by ambipolar diffusion and turbulent mixing. The 
lines b and c cannot escape, (d) The vertical ambipolar diffusion velocities act in opposite 
directions on a line of force to shorten its horizontal extent in the disk. They cannot remove 
it from the disk. 

Fig. 2 The variation of the magnetic field strength of a line of force for different initial conditions. 
The field strength at first grows by stretching, and later decreases because of ambipolar 
diffusion. The saturated state of the lines is the same, independent of its initial field strength. 
The initial values of B in microgauss were, .01, .03, .1, .3, 1.0, and 3. 

Fig. 3 The interstellar medium consists of dense clouds shown as spheres. A magnetic field line 
is shown which threads several clouds. The bowing up of the line in between the clouds is 
caused by cosmic ray and magnetic pressure in the intercloud medium. These pressures lead 
to a upward force on the neutrals in the cloud. The cloud is held down by gravitational 
force on the neutrals. The action of these two forces produces a motion of the ions through 
the cloud, ambipolar diffusion. The right hand end of the line is connected to the external 
cosmos. When the end of the line escapes from the last cloud by ambipolar diffusion, this 
part of the line will be expelled into the external world, and the part of the line in the disk 
will become shorter. 

Fig. 4 The lines of force for a nonuniform magnetic field after compression into the disk. To 
preserve zero divergence the lines of force must fan out in the direction of weaker field. It 
is seen that the region of weaker field, where B r is negative, is actually more extensive than 
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the region where B r is positive. Thus, after the field everywhere reaches a saturated field 
strength this weaker region will lead to a region of larger area, and the final mean magnetic 
field will be dominated by it. (Because the saturated value is independent of the initial field 
strength, the weakness of the initial field will not be important, and only its more extensive 
area will be important.) 

Fig.5 The variation of the magnetic field strength with z for different times as a result of the one 
dimensional calculation. At first the strength grows by stretching, and afterwards decreases 
by ambipolar diffusion. The value is negative because the initial radial field was positive, and 
the differential rotation sweeps it backwards. The magnetic field is plotted in dimensionless 
units. The dashed curves correspond to the times during which the field strength is 
increasing while the solid lines correspond to the decreasing phase. The convergence of the 
curves for the larger times is evident. 

Fig. 6 The variation of the magnetic field with time for initial conditions B r = B, L cos 9 for 9 
equal to 180 degrees in 15 degree increments, where Bi = .01 in dimensionless units. 
These curves are derived from the one dimensional calculation. The convergence of the 
curves at large t parallels the same convergence in figure 2, and the curves for equal initial 
conditions agree quantitatively. (The values of the initial conditions in microgauss are 
0.0285 x cos 9 + 0.000285 = 0.0288,0.0250,0.0205,0.0146,0.0077, etc. ) 

Fig. 7 The results of the two dimensional numerical simulation, for the case of an initial uniform 
magnetic field, Bq is shown as function of u at z = for different times. That the integral of 
Bq is zero is evident. This field will not lead to any Faraday rotation. The labeling of the 
curves is the same as in figure 5. 

Fig. 8 The results of the two dimensional numerical simulation, for the case of an initial non 
uniform magnetic field. Bq is shown as function of u at z = for different times. The regions 
of positive Bq are broader than those of negative Bq, so that the integral does not vanish, 
and this field will lead to a non zero Faraday rotation The labeling of the curves is the same 
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as in figure 5. 

Fig. 9 The results of the two dimensional numerical simulation, for the case of an initial non 
uniform magnetic field. Be is shown as function of z at u = for different times. The 
labeling of the curves is the same as in figure 5. The similarity to figure 5 is noteworthy. 

Fig. 10 The results of the two dimensional numerical simulation, for the both cases. Be is 

shown as function of t at u = 0, z = 0. The behavior is quite similar to figure 2 for the one 
dimensional simulation. 
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